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In conical refraction, when a focused Gaussian beam passes along one of the optic axes of a biaxial 
crystal it is transformed into a pair of concentric bright rings at the focal plane. We demonstrate both 
theoretically and experimentally that this transformation is hardly affected by partially blocking the 
Gaussian input beam with an obstacle. We analyze the influence of the size of the obstruction both 
on the transverse intensity pattern of the beam and on its state of polarization, which is shown to 
be very robust. 

PACS numbers: 42.25.Ja,42.25.Fx,78.20.Fm,42.25.Gy 


I. INTRODUCTION 

Gaussian beams are the most well known solution of 
the paraxial wave equation. They are form-invariant 
beams, i.e., the form of their transverse intensity pattern 
does not change upon propagation, apart from a scal¬ 
ing factor. Durnin et al. [1] reported another solution of 
the paraxial wave equation, the Bessel beams, which are 
completely invariant upon propagation. In other words, 
both the transverse intensity profile and scale of Bessel 
beams remain unchanged as it propagates, i.e., Bessel 
beams are diffraction-free beams. One of the main fea¬ 
tures of Bessel beams is that they self-reconstruct after 
an obstacle, being this effect known as self-healing. 

Recently, there has been a great interest in the study of 
the self-healing effect appearing in Bessel beams [2, 3] and 
other diffraction-free beams including Airy beams [4, 5] 
and Pearcey beams [6], or other exotic beams such as 
helico-conical beams [7] as well as Mathieu and Webber 
beams [8]. The major advantage of self-healing beams is 
that they can be used through turbulent media [5] and 
that they are ideal candidates for particle manipulation 
at different planes [9, 10] and in microscopy [11]. 

In conical refraction (CR) [12-14], when a circularly 
polarized focused Gaussian beam propagates through a 
biaxial crystal parallel to one of the optic axes it ap¬ 
pears transformed into a pair of concentric bright rings 
split by a null intensity (Poggendorff) ring at the focal 
plane, see Fig. 1(a). At this plane, the state of polar¬ 
ization of the rings is linear at every point, with the az¬ 
imuth rotating continuously so that every two diametri¬ 
cally opposite points have orthogonal polarizations. This 
non-factorization of the state of polarization makes CR 
beams a very particular type of vector beams [15]. Along 
the axial direction, the transverse intensity pattern of the 
CR beam becomes more involved and far enough from 
the focal plane an axial (Raman) spot is formed. The 
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above named features of CR beams are valid under the 
assumption that the CR geometrical ring radius Rq is 
much bigger than the waist radius of the focused input 
beam wq, i.e., assuming that po = Ro/^o 1- Fen 
smaller values of po, different light structures are formed 
[16-18]. 

CR has been reported as an efficient tool to generate 
Bessel beams [19-21]. The relation of CR with Bessel 
beams suggests that even if the input Gaussian beam 
is blocked partially, the CR beam may only be slightly 
affected. In this work we address the question of the 
CR healing of a Gaussian input beam in the presence 
of an obstruction. We investigate the reconstruction of 
the transverse intensity pattern of the CR beams at dif¬ 
ferent propagation distances and also of their state of 
polarization. We show that, even for relatively large ob¬ 
stacles, the CR beams keep their annular shape, state of 
polarization and dark optical singularities. The work is 
organized as follows: the theory of CR is presented in 
Section II. In Section III we describe our experimental 
set-up and report healing of conically refracted Gaussian 
beams after an obstruction for p 0 1. In Section IV we 
compare the reconstructing behavior of the CR beam for 
Po = 0.92 and po 1 and discuss the differences between 
the two cases. Finally, our main results are summarized 
in Section V. 



FIG. 1. (Color online) Transverse intensity patterns and state 
of polarization at the focal plane of CR beams numerically cal¬ 
culated from Eqs. (l)-(8) for a left handed circularly polarized 
Gaussian input beam for (a) po = Ro/wo = 10, (b) po = 1.5 
and (c) po as 0.92. Ro is the geometrical CR ring radius and 
wo the waist radius of the focused Gaussian input beam. 
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II. PARAXIAL THEORY OF CONICAL 
REFRACTION 

In what follows, we normalize our coordinates com¬ 
ponents to the beam waist wo and Rayleigh range zr , 
i.e., x x/wo , y — >• y/wo and 2 — >► zJzr. Additionally, 
note that when a light beam propagates through a biax¬ 
ial crystal (with three different refractive indices n x , n y 
and n z ) along one of the optic axis, there is a degeneracy 
on the refractive indices, i.e., n x = n z = n y = n [13]. 
Therefore, in the forthcoming equations only n needs to 
be taken into account (for the crystal used in the exper¬ 
iments described below n ~ 2 [14]). 

Let’s consider an input beam with transverse electric 
field 


E(r) — E x g x H - E y Q y , (l) 


where e x and e y are the unit vectors in Cartesian coor¬ 
dinates, and corresponding Fourier transform 


A(k) — Ac(k)e x + A y (k)e y7 

(2) 

OO 

A x (k) = JJ E x(*) e lkr dxdy, 

— OO 

(3) 

OO 

A/( k ) = J ^2 JJB y (r)e~ ik r dxdy. 

— OO 

(4) 

It can be shown that for a low birefringent biaxial crys¬ 
tal the beam evolution behind the biaxial crystal can be 
described as follows [25]: 

E x = B\' X {r, po) +B 0 ,y(r,p 0 ), 

E y = B 0 , x (r, po) + B 1}V (r, -p 0 ). 

(5) 

(6) 

with two main integrals 



So,a(r,Po) 

Bi,a{r,po) 


OO 

JJ e -l(k ' r_ ^ fe 2 )) h s in(po&) A a (k)dk x dky, 

— OO 
OO 

jJ^-2 JJ e -l(k ' r- ^ fe2) ^cos (p 0 k) + i ^ sin (p 0 h)j A a (k)dk x dk y , 

— OO 


(7) 

(8) 


where a = x,y. 

In this work, we consider a Gaussian input beam with 
electric field Ei n oc exp (— (x 2 + y 2 )) ei n , where ei n is 
the state of polarization of the beam. Fig. 1 shows the 
transverse intensity patterns and state of polarization at 
the focal plane of CR beams numerically calculated from 
Eqs. (l)-(6) for a left handed circularly polarized Gaus¬ 
sian input (eo = *4=(l,i)) beam for (a) po = 10, (b) 
po = 1.5 and (c) po = 0.92. As it can be observed, both 
the transverse intensity pattern and the polarization of 
the beams depend on po, as it has been shown in detail 
in [15]. For po = 10 1, all the features of the CR beam 

described above can be observed, namely: the transverse 
intensity pattern is formed by two bright rings split by 
the Poggendorff dark ring and every two opposite points 
are orthogonally linearly polarized. Note that under the 
condition po 1, the properties of the CR beam are un¬ 
affected by slightly different values of po- For Po = 1.5 the 
inner ring collapses into a point and the region separating 
the central spot and the outer ring possesses low but non¬ 
null intensity. For po = 0.92 the pattern is doughnut-like, 
being the central dark spot an exact null-intensity point, 
as reported in [17]. In both cases the state of polariza¬ 
tion is non-homogeneously distributed, including regions 
with elliptical states of different azimuth and ellipticity. 
Nonetheless, at the edge of the ring the polarization dis¬ 
tribution is linear with azimuth rotating an angle cj> for 


a full turn along the ring. The dependence of the CR 
transverse intensity pattern with po has been analyzed 
in detail in [16]. 


III. HEALING OF CONICALLY REFRACTED 
GAUSSIAN BEAMS FOR p 0 > 1 

In this Section we analyze the reconstruction of con¬ 
ically refracted Gaussian beams under the condition 
Ro wo when an obstruction of closing angle p blocks 
an azimuthal part of the input beam. We consider first 
the approximation Ro wo since it is the commonly 
used configuration in most experimental arrangements. 
Fig. 2 shows our experimental set-up. We obtain a cir¬ 
cularly polarized Gaussian beam at 640 nm from a diode 
laser coupled to a monomode fiber by utilizing, a lin¬ 
ear polarizer (LP) and a quarter wave-plate (QWP). As 
an obstruction, we use amplitude angular masks (AAM). 
They block an azimuthal sector of angle p of the Gaus¬ 
sian beam at the exit of the collimator. The AAM are 
made by printing the desired 2D pattern over a trans¬ 
parent sheet of plastic. The beam, whose waist radius 
wo can be adjusted by means of the collimator, passes 
through a biaxial crystal and parallel to one of the op¬ 
tic axes and a CCD camera combined with an imaging 
lens (IL) records the transverse intensity pattern of the 
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FIG. 2. (Color online) Experimental set-up. A Gaussian 
beam is obtained from a diode laser coupled to a monomode 
fiber. A collimator (C) is used to change the beam waist ra¬ 
dius wo of the Gaussian beam. A linear polarizer (LP) and a 
quarter waveplate (QWP) are used to fix the state of polar¬ 
ization of the Gaussian beam to circular. Azimuthal angular 
amplitude masks (AAM) characterized by their closing angle 
fi are used to block a sector of the Gaussian beam. The ob¬ 
tained beam passes through the biaxial crystal (BC) along one 
of its optic axes and a CCD camera and an imaging lens (IL) 
record the transverse intensity pattern at different planes. 


CR beam at different planes. As a biaxial crystal we 
use a commercially available (CROptics) KGd(WC> 4)2 of 
length l = 28 mm, conicity a = 16.9 mrad and, there¬ 
fore, Ro = la ~ 475 fi m. The polished entrance surfaces 
of the crystal (cross-section 6x4 mm 2 ) have parallelism 
with less than 10 arc sec of misalignment, and they are 
perpendicular to one of the two optic crystal axes within 
1.5 mrad misalignment angle for the central wavelength 
633 nm. The crystal is mounted over a micro-positioner. 
Fine alignment is performed by modifying independently 
the 0 and fi angles of the micro-positioner and observ¬ 
ing intensity pattern at the focal plane with the CCD 
camera. A cylindrically symmetric pattern is obtained 
when the input beam of circular polarization passes ex¬ 
actly along one of the optic axis of the BC. By reducing 
the focused beam waist down to wq = 44 gm, we have 
obtained a po parameter up to 10.75. 

Fig. 3 shows both the experimental (top row) and nu¬ 
merically calculated (bottom row) transverse intensity 
pattern at the focal plane for obstructions of (a,f) 0 = 0°, 
(b,g) 45°, (c,h) 90°, (d,i) 135° and (e,j) 180°. Insets rep¬ 
resent the obstructed input beam just behind the AMM. 
When a relatively small obstruction angle is considered 
(fi = 45°), the transverse intensity pattern is almost un¬ 
affected as compared with the case without obstruction. 
In this case, the Poggendorff dark ring and the two ring¬ 
like structures are clearly visible. The rings are mirror- 
symmetric with respect to the horizontal axis but they 
have a maximum at their top and bottom regions. As 
< f> increases the CR rings become asymmetric and the 
outer ring breaks into two boomerang-like lobes such that 
no complete Poggendorff dark ring is appreciable. For 
< j) = 180°, i.e., when the AAM blocks half of the input 
beam, the CR transverse intensity pattern is formed by a 
wide single ring with two dark singularities in the upper 


and top regions of the ring. 




FIG. 3. (Color online) Transverse intensity patterns at the 
focal plane (z — 0) obtained with an obstruction of closing 
angle (a,f) </> = 0°, (b,g) fi = 45°, (c,h) fi = 90°, (d,i) fi = 
135° and (e,j) fi = 180° placed before the biaxial crystal. 
(a,f) Intensity pattern obtained in the absence of obstruction. 
First row: experimental results (/9 q XP — 10.75). Second row: 
numerical calculations obtained from Eqs. (l)-(6) (po h = 10). 
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FIG. 4. (Color online) Numerically calculated 2D density 
plots of the azimuth e (first row) and the ellipticity fi (second 
row) of conically refracted beams at the focal plane for po = 10 
when the input Gaussian is blocked by azimuthal obstructions 
of (b,g) fi = 45°, (c,h) fi — 90°, (d,i) fi = 135° and (e,j) 
< j> = 180° is placed before the biaxial crystal. (a,f) 2D density 
plots of e and fi in the absence of the obstruction. 

The healing mechanism provided by the biaxial crys¬ 
tal to reconstruct the CR beam can be understood in 
terms of the wave-vector splitting within the crystal. 
Every plane-wave is described by a certain wave-vector 
k = k|| +kx, with k^ = |&_l|(cos fik, sinfik). The biax¬ 
ial crystal splits every plane wave into two new plane- 
waves. At the focal plane, these two plane-waves are 
refracted at positions on the ring characterized by their 
azimuthal angle p = fik and p = fik + tt [22]. As a 
consequence, when one azimuthal sector of the Gaussian 
beam is blocked, the azimuthally opposite sector partially 
compensates the absence of the blocked sector. For this 
reason, even when half of the input beam is blocked, a 
ring-like structure can be formed after passing through 
the biaxial crystal. 

This mechanism explains why a single bright ring with¬ 
out Poggendorff splitting is obtained when half of the in¬ 
put beam is blocked with the AAM, i.e., when fi = 180°. 
The two bright rings with Poggendorff splitting appear 
as an interference of plane waves going to a particular 
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azimuthal point of the CR ring from opposite sectors of 
the input beam. In contrast, if there are no other waves 
coming to the corresponding opposite points of the CR 
ring pattern at the focal plane, there is no interference 
and a only a single ring is observed. In Figs. 3(b)-(e) 
it is clearly visible an increase of the azimuthal sector of 
single bright ring and shrink of the double bright rings’ 
domain with the AMM closing angle 0. Note that the 
azimuthal sector occupied by the double bright rings is 
larger than the angular sector of the AMM, 0, because of 
the diffraction of the input beam at the edges of AMM 
dark sector. 

Now we turn to analyze the healing of the state of po¬ 
larization of the reconstructed CR beams. The standard 
tool to analyze the state of polarization of a light beam 
is the Stokes vector: S = (So, Si, $2, S 3 ). ^or an e l ec tric 
field E = (. E x ,E y ) with intensity I the Stokes parameters 
read [26]: 


So = I= \E x \ 2 + \Ey\ 2 , 

(9) 

Si = J 0 O — J9OO = \E X \ — \Ey\ , 

(10) 

S 2 = /45° - /1350 = 2Re [E*E y ], 

(11) 

S 3 = Ir~Il= 21m [E*E y ], 

(12) 


where (<F = 0°,45°,90°, 135°) indicates the intensity 
of linearly polarized light with azimuth <F, and Ir and 
II indicate the intensity of right- and left-handed circu¬ 
larly polarized light, respectively. In what follows, we use 
equations normalized to E 2 , i.e., we consider I = E 2 = 1. 
These definitions show that So account for the intensity 
of the light beam, Si measures the amount of light which 
is linearly polarized (LP) in the vertical/horizontal ba¬ 
sis, S 2 does the same but with the diagonal basis and 
S 3 relates the state of polarization in the right- and left- 
circular ly polarized (CP) basis. The following equations 
show how the Stokes parameters account for the azimuth 
e and ellipticity /3 of the polarization ellipse [26] : 


1 ! 

e = - arctan 

2 ' 

kSiJ 

(13) 

/3 = ^ arctan j 

( s 3 \ 

(14) 

W s ‘i + s‘i) 


The values of e and j3 at the focal plane numerically cal¬ 
culated for obstructions of (a,f) 0 = 0°, (b,g) 45°, (c,h) 
90°, (d,i) 135° and (e,j) 180° are presented in Fig. 4. 
Insets represent the obstructed input beam. Ideally, as 
po 00 , [3 0. Since in our numerical simulations 

we consider po = 10, the state of polarization of the CR 
rings is slightly elliptical rather than purely linear. As 
it can be appreciated, in general, the polarization struc¬ 
ture of the CR beam is maintained for all the values of </>, 
i.e., every two diametrically opposite points of the light 
structure are orthogonally polarized. 

In Fig. 5 we show the evolution of the transverse in¬ 
tensity pattern along the axial direction for a blocking 
sector of <j) = 45°. Top row shows the transverse inten¬ 
sity patterns in the absence of the blocking mask, while 


middle and bottom rows are the experimental and nu¬ 
merically calculated transverse intensity patterns for the 
obstructed Gaussian beam. Near the focal plane (z = 0) 
the transverse intensity pattern resembles the pattern ob¬ 
tained without obstruction. In contrast, far enough of 
the focal plane it can be appreciated a perturbation of 
the CR transverse intensity pattern that resembles the 
considered obstruction. 


IV. RECONSTRUCTION PROPERTIES OF 
CONICALLY REFRACTED GAUSSIAN BEAMS 

FOR po ~ 1 

As discussed in Section II, the CR beam depends 
strongly on the value of the control parameter po. For 
po = 0.92 the transverse intensity pattern at the focal 
plane forms a doughnut-like light structure with a null 
intensity point at the beam center [17], see Fig. 6(a). 
Along the axial direction, the intensity at the beam cen¬ 
ter is no longer zero and the beam forms an optical bottle. 
This value of po is particularly interesting because the 
polarization distribution of the light ring has points with 
different /3 and e and the beam forms a Poincare beam, 
i.e., a beam possessing points with all the polarization 
states of the Poincare sphere [27]. In what follows we 
discuss the reconstruction of a conically refracted Gaus¬ 
sian beam for po =0.92 after an obstruction, analogously 
to what has been performed above for po 1. For the 
incoming experiments we have used the same set-up as 
in Fig. 2 but with a 2.3 mm long KGd(W 04)2 (therefore 
Ro = 39 /im) crystal and a waist radius of wo = 44 /im. 

Fig. 6 shows both the experimental (top row) and nu¬ 
merically calculated (bottom row) transverse intensity 
pattern at the focal plane for obstructions of (a,f) 0 = 0°, 
(b,g) 45°, (c,h) 90°, (d,i) 135° and (e,j) 180°. Insets 
represent the obstructed input beam for p 0 = 0.92. In 
contrast to the case of po 1 , when a relatively small 
obstruction angle is considered (0 = 45°), the transverse 
intensity pattern is substantively different with respect 
to the case with no obstruction. In this case, a maxi¬ 
mum of intensity appears at the bottom part of the ring 
and, therefore, the intensity pattern is mirror symmetric 
with respect to the vertical axis. As </> increases, the in¬ 
tensity in the bottom part of the light structure becomes 
stronger than in the top part. For all the values of 0, 
an intensity minimum can be observed but its position 
moves in the vertical direction as </> increases. However, 
note that even when half of the beam is blocked by the 
obstruction, the dark singularity is preserved. A similar 
behavior has been reported for a linearly polarized Bessel 
beam [3]. 

With respect to the state of polarization of the CR 
beams at the focal plane, see Fig. 7, we have observed 
that there is a tendency to preserve the polarization 
structure of the CR beam without obstruction: the CR 
beam is elliptically polarized and the ellipticity changes 
radially similarly to the intensity pattern. At the edges of 
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FIG. 5. (Color online) Experimental (second row) and numerically calculated with Eqs. (l)-(6) (third row) transverse intensity 
patterns along the axial direction z for an obstruction with 0 = 45 for po = 10. The first row shows the numerically calculated 
transverse intensity patterns obtained in the absence of the obstruction, z is measured in units of the Rayleigh range, which 
for the Gaussian beam used in our experiments (wo = 44 pm) is zr — 9.5 mm. 



FIG. 6. (Color online) Transverse intensity patterns at the 
focal plane (z = 0) obtained when an obstruction of angle 
(b,g) 4> = 45°, (c,h) 4> = 90°, (d,i) 0 = 135° and (e,j) 4> = 180° 
is placed before the biaxial crystal. (a,f) Transverse intensity 
pattern in the absence of obstruction. First row: experimental 
results (Pq XP — 1-04). Second row: numerical calculations 
obtained from Eqs. (l)-(6) (po h = 0.92). 


FIG. 7. (Color online) 2D density plots of the azimuth e (first 
row) and the ellipticity /3 (second row) of conically refracted 
beams at the focal plane for po = 0.92 when the input Gaus¬ 
sian is blocked by azimuthal obstructions of (b,g) 4> = 45°, 
(c,h) <f> = 90°, (d,i) <f> = 135° and (e,j) 0 = 180° is placed 
before the biaxial crystal. (a,f) 2D density plots of e in the 
absence of obstruction. 


the beam, /? —>► 0 and the characteristic CR polarization 
distribution is recovered. However, since the transverse 
intensity pattern is very affected by the presence of an 
obstruction, these features of the state of polarization are 
lost for large enough values of <fi. In particular, (3 losses 
its doughnut-like shape, while e is kept quite stable up 
to (j) = 135°. 

The evolution of the transverse intensity pattern along 
the axial direction for a blocking sector of 0 = 45° is 
shown in Fig. 8. Top row shows the transverse inten¬ 
sity patterns in the absence of the blocking mask, while 
middle and bottom row are the experimental and numer¬ 
ically calculated transverse intensity patterns for the ob¬ 
structed Gaussian beam. Near the focal plane (z = 0) the 


transverse intensity pattern differs from the one obtained 
without obstruction. Additionally, out of the focal plane 
there is no reconstruction of the transverse intensity pat¬ 
tern. Therefore, for po = 0.92, the beam reconstruction 
process is not found at the focal plane nor away from it. 


V. CONCLUSIONS 

We have analyzed in detail the transformation of Gaus¬ 
sian beams partially obstructed when they propagate 
through a biaxial crystal and parallel to one of the op¬ 
tic axes, i.e., under conditions of conical refraction (CR). 
We have shown that, at the focal plane, the CR beams 
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FIG. 8. (Color online) Experimental (second row) and numerically calculated with Eqs. (l)-(6) (third row) transverse intensity 
patterns along the axial direction z for an obstruction with <f> — 45 for po = 0.92. The first row shows the numerically calculated 
transverse intensity patterns obtained in the absence of the obstruction, z is measured in units of the Rayleigh range, which 
for the Gaussian beam used in our experiments (wo = 44 pm) is zr — 9.5 mm. 


for po 1 preserve the annular shape even when half 
of the beam is blocked. However, we have found that 
the dark annular singularity known as Poggendorff dark 
ring only remains for small perturbations of the input 
beam. Out of the focal plane we have obtained that the 
obstruction affects the beam evolution, being its effect 
more appreciable the further one moves along the axial 
direction. Additionally, the reconstruction of the state 
of polarization of the CR beam has also been investi¬ 


gated. For po 1, we have found that the polarization 
distribution of the CR rings is very stable against large 
perturbations. 

We have carried out analogous investigations for po ~ 
1. In this case the transverse light pattern is more af¬ 
fected by the presence of the obstruction than in the case 
of po 1. Regarding the reconstruction of the state of 
polarization we have found that only the azimuth of the 
polarization is relatively robust when large obstructions 
affect the input Gaussian beam. 
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